Fault-Tolerant Quantum Dynamical Decoupling 
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Dynamical decoupling pulse sequences have been used to extend coherence times in quantum 
systems ever since the discovery of the spin-echo effect. Here we introduce a method of recursively 
concatenated dynamical decoupling pulses, designed to overcome both decoherence and operational 
errors. This is important for coherent control of quantum systems such as quantum computers. For 
bounded-strength, non-Markovian environments, such as for the spin-bath that arises in electron- 
and nuclear-spin based solid-state quantum computer proposals, we show that it is strictly advanta- 
geous to use concatenated, as opposed to standard periodic dynamical decoupling pulse sequences. 
Namely, the concatenated scheme is both fault-tolerant and super-polynomially more efficient, at 
equal cost. We derive a condition on the pulse noise level below which concatenated is guaranteed 
to reduce decoherence. 
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Introduction. — In spite of considerable recent progress, 
coherent control and quantum information processing 
(QIP) is still plagued by the problems associated with 
controllability of quantum systems under realistic condi- 
tions. The two main obstacles in any experimental re- 
alization of QIP arc (i) faulty controls, i.e., control pa- 
rameters which are limited in range and precision, and 
(ii) decoherence- errors due to inevitable system-bath in- 
teractions. Nuclear magnetic resonance (NMR) has been 
a particularly fertile arena for the development of many 
methods to overcome such problems, starting with the 
discovery of the spin-echo effect, and followed by meth- 
ods such as refocusing, and composite pulse sequences 
. Closely related to the spin-echo effect and refocusing 
is the method of dynamical decoupling (DD) pulses intro- 
duced into QIP in order to overcome decoherence-errors 
1^ 1^ . In standard DD one uses a periodic sequence of 
fast and strong symmetrizing pulses to reduce the unde- 
sired parts of the system-bath interaction Hamiltonian 
HsB, causing decoherence. Since DD requires no encod- 
ing overhead, no measurements, and no feedback, it is an 
economical alternative to the method of quantum error 
correcting codes (QECC) [e.g., 0, H, Hi and references 
therein] , in the non-Markovian regime [7| . 

Here we introduce concatenated DD (CDD) pulse se- 
quences, which have a recursive temporal structure. We 
show both numerically and analytically that CDD pulse 
sequences have two important advantages over standard, 
periodic DD (PDD): (i) Significant fault-tolerance to 
both random and systematic pulse-control errors (see 
Ref. [ll for a related study), (ii) CDD is significantly 
more efficient at decoupling than PDD, when compared 
at equal switching times and pulse numbers. These ad- 
vantages simplify the requirements of DD (fast-paced 
strong pulses) in general, and bring it closer to utility 
in QIP as a feedback-free error correction scheme. 

The noisy quantum control problem. — The problem of 



faulty controls and decoherence errors in the context of 
QIP, as well as other quantum control scenarios 01, can 
be formulated as follows. The total Hamiltonian H for 
the control-target system [S] coupled to a bath {B) may 
be decomposed as: H = Hs^Ib + Is'^Hb + Hsb, where 
/ is the identity operator. The component Hsb is respon- 
sible for decoherence in S. We focus here on the single 
qubit case, but the generalization to many qubits, with 
Hsb containing only single qubit couplings, is straight- 
forward. We shall interchangeably use X, Y, Z to denote 
the corresponding Pauli matrices ctq, and ctq to denote 
/. The system Hamiltonian is Hs — H™^ + Hp, where 
i/lj"* is the intrinsic part (self Hamiltonian), and Hp is 
an externally applied, time-dependent control Hamilto- 
nian. We denote all the uncontrollable time-independent 
parts of the total Hamiltonian by He, the "else" Hamil- 
tonian: He ■— H™^ + Hb + Hsb- We assume that all 
operators, except /, are traceless (traceful operators can 
always be absorbed into Hs (8> Iq, Is ^ Hb). We fur- 
ther assume that \\He\\ < oo [iflj. Note that this is 
a physically reasonable assumption, even in situations 
involving theoretically infinite-dimensional environments 
(such as the modes of an electromagnetic field), since in 
practice there is always an upper energy cutoff We 
consider "rectangular" pulses [piece- wise constant Hp{t)] 
for simplicity; pulse shaping can further improve our re- 
sults An ideal pulse is the unitary system-only op- 
erator P{6) = Texp[— i Hp{t)dt], where T denotes 
time-ordering and h — 1 units are used throughout. A 
non-ideal pulse, Up{S) = Tcxp[— i ^^{Hp{t) + Wp{t) + 
He{t)}dt\, includes two sources of errors: (i) Deviations 
Wp from the intended Hp. Such deviations can be ran- 
dom and/or systematic, generally operator- valued; (ii) 
The presence of He during the pulse. 

Periodic DD. — In standard DD one periodically ap- 
plies a pulse sequence comprised of ideal, zero-width 
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TT-pulses representing a "symmetrizing group" Q = 

{PiYiJo (-^0 = I), and their inverses. Let fi-o = 
T exTp[—i He{t)dt] denote the inter-pulse interval, 
i.e., free evolution period, of duration tq. The effec- 
tive Hamiltonian He^^ for the "symmetrized evolution" 
nUcT^ P/f ro-P; =: 6*1^1'^°-^=^' is given for a single cycle 
by the first-order Magnus expansion: H^p « H^s = 
I^EE'tT^ p}HP^ 2]. This result is the basis of an el- 
egant group-theoretic approach to DD, which aims to 
eliminate a given Hsb by appropriately choosing Q 2, ^ . 
The "universal decoupling" pulse sequence, constructed 
from C/uD {co, ci, (72, cs}, proposed in plays a cen- 
tral role: it eliminates arbitrary singie-qubit errors. For 
this sequence we have, after using Pauli-group identities 
{XY = Z and cyclic permutations), pi := 6*"^^^= ' = 

Y{l=oPhroP^ = ±roXfroZ±roXtroZ, wherC T, = 4to. 

The idea of dynamical symmetrization has been thor- 
oughly analyzed and appHed (see, e.g., 0, and ref- 
erences therein). However, higher-order Magnus terms 
can in fact not be ignored, as they produce cumulative 
decoupling errors. Moreover, standard FDD is unsuited 
for dealing with non-ideal pulses 8]. 

Concatenated DD. — Intuitively, one expects that a 
pulse sequence which corrects errors at different levels of 
resolution can prevent the buildup of errors that plagues 
FDD; this intuition is based on the analogy with spa- 
tially-concaXenated QECC (e.g., ji|). With this in mind 
we introduce CDD, which due to its temporal recursive 
structure is designed to overcome the problems associ- 
ated with FDD. 

Definition 1 A concatenated universal decoupling pulse 
sequence: p„+i := p„Xp„Zp„Xp„Z, where po = fro and 
n > 0. 

Several comments are in order: (i) pi is the "universal 
decoupling" mentioned above, but one may of course also 
concatenate other pulse sequences; (ii) One can interpret 
pi itself as a one-step concatenation: pi '.—pxYjixY , 
where px ■— fXfX (f-.^frg) and Fauli-group identities 
have been used, (iii) Any pair, in any order, of unequal 
Fauli TT-pulses can be used instead of X and Z, and fur- 
thermore a cyclic permutation in the definition of pi is 
permissible; (iv) The duration of each sequence is given 
by 7" ^ T-n '■— 4"to (after applying Pauli-group iden- 
tities); (v) The existence of a minimum pulse interval 
To and finite total experiment time T are practical con- 
straints. This sets a physical upper limit on the num- 
ber of possible concatenation levels rimax in a given ex- 
periment duration; (iv) Pulse sequences with a recursive 
structure have also appeared in the NMR literature (e.g., 
[l3|'). though not for the purpose of reducing decoherence 
on arbitrary input states. We next present numerical sim- 
ulations which compare CDD with FDD. 

Numerical Results for Spin-Bath Models. — For com- 
paring the performance of CDD vs FDD, we have cho- 



sen an important example of solid-state decoherence: a 
spin-bath environment This applies, e.g., to spec- 

tral diffusion of an electron-spin qubit due to exchange 
coupling with nuclear-spin impurities e.g., in semi- 
conductor quantum dots |16||, or donor atom nuclear 
spins in Si [17l |. Specifically, we have performed numer- 
ically exact simulations for a model of a single qubit 
coupled to a linear spin-chain via a Heisenberg Hamil- 
tonian: He = UJsCrf + UJB{Y.a=2 O + Ef>fc>l Jab^a ' ^b- 

The system spin-qubit is labelled 1; the second sum rep- 
resents the Heisenberg coupling of all spins to one an- 
other, with jab — J exp(— Adab), where A is a constant 
and dab is the distance between spins. Such exponen- 
tially decaying exchange interactions are typical of spin- 
coupled quantum dots The initial state is a ran- 
dom product state for the system qubit and the envi- 
ronment. The goal of DD in our setting is to minimize 
(the log of) the "lack of purity" of the system qubit, 
I = log]^o(l ~ Tr[p|]), where ps is the system density 
matrix obtained by tracing over the environment basis. 
At given CDD concatenation level n we also implement 
FDD by using the same minimum pulse interval tq as 
in CDD and the same total number of pulses iV < 4"; 
this ensures a fair comparison. In all our simulations we 
have set the total pulse sequence duration T = 1 , in units 
such that {ujsT,uibT, Xdjj+i) = (2, 1,0.7). Longer pulse 
sequences correspond to shorter pulse intervals tq. Note 
that we have chosen our system and bath spins to be sim- 
ilar species. Qualitatively, the number of bath spins K 
had no effect in the tested range 2 < K < 7, while quan- 
titatively, and as expected, decoherence rises with K. 
DD pulses were implemented by switching Hp ~ haf, 
a G {x,z}, on and off for a finite duration S > 0; note 
that n < nniax(T, (5). We define the pulse jitter Wp as 
an additive noise contribution to Hp. It is represented 
as Wp = r** • (Ji , with being a vector of random (uni- 
formly distributed) coefficients. We distinguish between 
systematic {Wp fixed throughout the pulse sequence, but 
different for each a) and random {Wp changing from 
pulse to pulse) errors. 

Our simulation results, shown in Figs. I1I3I compare 
CDD and FDD as a function of coupling strength, rel- 
ative jitter magnitude, and number of pulses. Fig. ^ 
left, compares CDD and FDD at fixed number of pulses. 
CDD outperforms FDD in the random jitter case with 
noise levels of up to almost 10%. Fig.^ right, shows the 
performance of CDD as a function of jitter magnitude 
and concatenation level: the improvement is systematic 
as a function of the number of pulses used. Figure[21con- 
trasts CDD and FDD in the jitter-free case, as a function 
of system- bath coupling j. As predicted in the analytical 
treatment below, CDD offers improvement compared to 
FDD in decoherence reduction over a wide range of j val- 
ues. Figure 13 compares CDD and FDD as a function of 
systematic jitter. Superior performance of CDD is par- 
ticularly apparent. These results establish the advantage 
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FIG. 1: Left: Performance of CDD (solid line, n = 4) vs 
PDD (dot-dashed line) as a function of random jitter fraction 
\wp\ ■- \\Wp\\ I lliypll, with pulse-width & = KT^T , coupling 
strength j — .2/T (T is the total evolution time), and number 
of bath spins K = 2, averaged over 90 jitter realizations. For 
comparison, the horizontal dashed line corresponds to free 
evolution. Right: CDD as a function of random jitter fraction 
Iwpl and concatenation level n. The vertical axis denotes 
I — logj^g(l — purity) here and in Figs. 2 and 3. 
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FIG. 2: CDD (left) and PDD (right), as a function of system- 
bath coupling j; pulse width 5 = 10~'^T, number of bath spins 
K = 5, and without jitter {Wp — 0). Note the Z-axis scale 
difference between CDD and PDD. 



of CDD over PDD in a model of significant practical in- 
terest, subject to a wide range of experimentally relevant 
errors. We now proceed to an analytical treatment. 

Imperfect decoupling. — Consider DD pulse sequences 
composed of ideal pulses. Let us partition He as He — 
Hji + h\, where Hj^ = Y ® By + Z ® B^, and h\ = 
X ® Bx + Hb- The super/subscripts ±, X and ||, X cor- 
respond to terms that anti-commute and commute with 
X ® Ib: respectively. Thus the effect of px =fXfX 
in PDD can be viewed as a projection of Hsb onto the 
component "parallel" to X, i.e., h\. For the Y pulses in 




FIG. 3: Left: CDD performance as a function of concatena- 
tion level and systematic jitter. The pulse width 5 = 10~*r, 
number of bath spins A" = 5, jT = 15.0, averaged over 7 — 80 
realizations (more realizations for higher n). Right: CDD 
(solid line) vs PDD (dot-dashed line) as a function of system- 
atic jitter for n = 5, 5 = 10~^T , if = 5, jro = 3.0, averaged 
over 14 — 80 realizations. The dashed line is pulse- free evo- 
lution. CDD performance is unaffected up to ~ 20% jitter 
level. 
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FIG. 4: Projections involved in DD. (a) Perfect cancelation 
in first-order Magnus case, (b) Extra rotation induced by 
higher-order Magnus terms, and efi'ect of concatenation. 



pi =^xYvxY we can similarly write H\ = Hy + Hy, 
where -L,Y (||,i^) denotes anti-commutation (commuta- 
tion) with Y, whence Hy = Hb- Then the role of the 
Y pulses is to project onto Hy, which eliminates 
Hsb altogether, i.e., transforms He = Hsb + Hb into 
a "pure-bath" operator Hb- This geometrical picture 
of two successive projections is illustrated in Fig. QJa). 
However, these projections are imperfect in practice due 
to second-order Magnus errors. Indeed, instead of a 
sequence such as fXfX, one has, after pulse Pi=x,Y, 
Ie,, exphzT(i?P - Ht)]eM-^<H} + H^)], where 
He — nf + H^, and we have accounted for the sign- 
flipping due to Pi- Using the BCH formula (e.g., [131 )j 
we approximate the total unitary evolution as lE,i = 
exp[-z(2T)i7eff,,; + O (Af)], where iJeff,, :=I?p. (T)[He] = 

^-^rHt/2^\\^^rHt/2^ whcrC Af 1= T^\\H^f\\H}\\, and it 

is assumed that, since \\He\\ < oo, one can pick r such 
that Ai ^ 1. The mapping 'Dp.{T)[He] clearly has a geo- 
metric interpretation as a projection that eliminates H^^, 
followed by a rotation generated by H^^- This rotation 

produces extra system-bath terms besides Hf , hence im- 
perfect DD- This is illustrated in Fig. EJb) : in the first- 

(2) 

order Magnus approximation the transition from Hg^ 

0. But 

(3) 

in the presence of second-order Magnus errors iJgg 0. 
The difference between CDD and PDD is precisely in the 

(3) 

manner in which this error is handled: in PDD the H^"^ 



to i/gg suffices to eliminate -ffsB, i-e., H^ 



error accumulates over time since the same procedure is 
simply repeated periodically. However, in CDD the pro- 
cess of projection-|-rotation is continued at every level of 
concatenation, as suggested in Fig.^Jb) (red arrow above 
H^g^2)- In CDD, i/^g' is shrunk with increasing m, m a 
manner we next quantify. 

Convergence of CDD in the limit of zero-width 
pulses. — Decoupling induces a mapping on the com- 
ponents of He- For a single qubit, writing He = 



^ a—x ,y ,z 



Bn. we have He 



pi 



H. 



(1) 



B^'a^ , where a second-order Magnus expansion yields: 
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= Bo, Bi'^ = iTo[Bo,B,], b\^^ = iTo\{[Bo,By] - 
i{B^,B^}), Bi^^ ^ 0. Let us define /3 := ||Bo|| 
and J :— max{\\Bx\\ ,\\By\\ ,\\Bz\\), where we assume 
J < (3 < oo Comparing with the model we 

have used numerically, J = 0{\) and /3 ~ 0{ujb)- 
It is possible to show that a concatenated pulse se- 
quence p„ can still be consistently described by a second- 
order Magnus expansion at all levels of concatenation, 
provided the (sufRcient) condition t„/3 <g; 1 is sat- 
isfied [i^l- We can then derive the recursive map- 
ping relations for iji""^-* He^^ = J^a'^a ® Ba^^: 
5'"^°) = Bo, Bi"^') = (jr„_i)[Bo,Bi"-'^], 4""'^ = 
i(iT„_i)[i?o, =0. The propagator corre- 

sponding to the whole sequence is ex\){~iTnHe^^), which 
in the limit of ideal performance reduces to the identity 

in) 

operator. These results for Ba allow us to study the 
convergence of CDD, and bound the success of the DD 
procedure, as measured in terms of the fidelity (state 
overlap between the ideal and the decoupled evolution). 
This fidelity is given by Q 

/„ « l-||r„ij(")||2 « l-(r„/i("))2 =: l-(<i>cDD)', (1) 

where H is the system-traceless part of H, and h^^^ :— 
max{ I \Bi"-'> 11,1 1 1 }. We find that 

$CDD < iPT/N^/^riJT), (2) 

where T = Ntq < Tn = 4"to is the total sequence 
duration, comprised of N pulse intervals. In contrast, 
$PDD = yields 

$PDD = 2(/3ro)(jr) = 2(/3r/iV)(JT). (3) 

Note that for = 4, $cdd = 'I'pdd as expected. There 
is a physical upper limit to the number of concatena- 
tion levels, imposed by the condition /3r„ ^ 1. Us- 
ing this condition in the form /3 = c/T, where c is 
some small constant (such as 0.1), and fixing the value 
of P, we can back out an upper concatenation level 
f^max — — log4 ^ ; inserting this into Eq. (0) we have 

*I>CDD < (c/3ro)^^ ^ (JT). We can now compare 
the CDD and PDD bounds in term of the final fidelity: 

1 - /cDD (c/3ro)-'°S4^ pro^o 

1 - /pDD - 4(/?To)2 ^ ^ > 

This key result shows that CDD converges super- 
polynomially faster to zero in terms of the (physically 
relevant) parameter /^tq, at fixed pulse sequence duration. 
However, it is important to emphasize that our bound on 
'I'cDD is unlikely to be very tight, since we have been very 
conservative in our estimates (e.g., in applying norm in- 
equalities and estimating convergence domains). Indeed, 
in our simulations (above) /3t„ « 2, which is beyond our 
conservatively obtained convergence domain. 



Finite width pulses. — We now briefly consider 
the more realistic scenario of rectangular pulses 
r exp[-i Jo{Hp{t) + He{t)}dt] of width S < tq. In 
this case we can derive a modified form of the condition 
/3r„ ^ 1, required for consistency (of using a second- 
order Magnus expansion at all levels of concatenation) 
113: 

cr„/3-hd— <Cl, (5) 

Tn 

where c, d ~ 1 are pulse sequence-specific numerical fac- 
tors. The consistency requirement (jSJ validates the anal- 
ysis of convergence of CDD for S ^ 0, and we can repro- 
duce the advantage of CDD over PDD for 5 = [mani- 
fest in Eq. Q]. As expected Eq. © imposes a more de- 
manding condition on the total duration t„, at fixed bath 
strength f3. While Eq. ^ cannot be called a threshold 
condition (in analogy to the threshold in QEC), since it 
depends on the total sequence duration, it does provide 
a useful sufficient condition for convergence of a finite 
pulse-width CDD sequence, and introduces the concept 
of error per gate which is fundamental in QEC. 

Conclusions and outlook. — We have shown that con- 
catenated DD pulses offer superior performance to stan- 
dard, periodic DD, over a range of experimentally rele- 
vant parameters, such as system-bath coupling strength, 
and random as well as systematic control errors. Here 
we have addressed the preservation of arbitrary quan- 
tum states. Quantum computation can in principle be 
performed, using CDD, over encoded qubits by choosing 
the DD pulses as the generators of a stabilizer QECC, 
and the quantum logic operations as the corresponding 
normalizer [21I . Another intriguing possibility is to 
combine CDD and high-order composite pulse methods 
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